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Tho  program  of  research  offered  in  the  work  proposul  for  P.Y.  11)70 
(December  1,  1969  -  November  30,  1970)  has  been  essentially  fulfilled  as 


follows: 

1)  We  have  established  a  suitable  Green’s  function  model  to  explain 
the  generation  of  long  period,  large  amplitude  gravity  modes  by 
thermonuclear  and  nuclear  tests  in  the  atmosphere.  Thus,  it  is 
shown  in  a  paper  in  press  by  Tolstoy  and  Lau^,  that  the  impor¬ 
tant  mechanism  here  must  be  the  rising  sphere  of  hot  gas.  It  is 
not  difficult  to  demonstrate  that  -  for  large  yields,  at  least  - 
this  mechanism  should  be  orders  of  magnitude  more  efficient  in 
generating  long  period  internal  gravity  wave  modes  than  the  ex¬ 
plosion  phase  proper.  It  is  our  belief  now  that  this  is  the 
source  of  observed  long  wavelength  gravity  waves;  although  we 
have,  in  the  attached  report,  merely  illustrated  the  principle 
of  the  mechanism,  without  extending  the  results  to  a  realistic 
model  of  the  atmosphere,  there  appears  to  be  little  doubt  that 
one  may  in  this  manner  account  for  the  large  amplitudes  of  ver¬ 
tical  motion  observed  for  such  disturbances  traveling  in  the. 
ionosphere . 

2)  A  further  product  of  the  undersigned's  investigations  has  been 

embodied  in  another  article  entitled  "Damping  of  very  long  gravity 

(2) 

waves  in  tho  atmosphere"'  y:  In  this  work,  tho  question  of  inter¬ 
nal  gravity  wavo  damping  has  boon  examined,  based  upon  the  assump- 
tion  that  most  of  the  attenuation  occurs  in  a  sort  of  atmospheric 
surface  layer  in  which  the  mean  free  paths  of  the  neutral  gas 
molecules  are  comparable  to  the  wavelengths  involved.  This  layer 
may  be  viewed  as  a  transition  between  two  regimes:  the  effective 
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vacuum  above  and  the  region  below,  for  which  the  concepts  of  con¬ 
tinuum  mechanics  apply.  Simple  kinetic  theory  arguments  predict 
then  the  reflexion  of  long  internal  waves  from  this  layer,  the 
reflexion  coefficient  being  about  e"a^  where  a  is  a  constant  of 

the  order  of  it,  y  is  the  vertical  wavenumber  and  d  the  local  scale 

(31 

height.  This  turns  out  to  be  in  accord  with  Yanowitch's  result 
for  the  reflexion  of  internal  gravity  waves  by  a  region  of  increas 
ing  viscosity  in  an  isothermal  atmosphere,  a  result  obtained  from 
asymptotic  arguments  applied  to  the  equations  of  continuum  mechan¬ 
ics.  Thus  it  appears  possible  to  account  for  a  certain  number  of 
observations  of  long  internal  gravity  wave  observations  traveling 
to  great  distances  (e.g.,  round  the  world  paths)  as  waveguide 
phenomena,  and  to  justify  the  application  of  free  surface  bound¬ 
ary  conditions  to  the  atmosphere. ^4* 5 

3)  The  problem  of  deducing  explosion  diagnostics  from  these  distur¬ 
bances  has  also  been  confronted;  a  discussion  of  this  problem 
has  been  embodied  in  a  previously  circulated  report where  it 
was  concluded  that  the  chief  practical  obstacle  to  the  use  of 
observed  disturbances  in  this  manner  is  a  lack  of  sufficient 
understanding  of  the  detailed  properties  of  the  upper  atmosphere. 

A  method  of  remedying  this  situation  hi'S  i'een  discussed  in  that 
(71  • 

report  ,  and  was  proposed  in  our  F.Y.  1971  work  proposal. 

4)  Of  immediate  relevance  to  the  last  item  have  been  studies  carried 

(81 

out  jointly  by  this  '.riter  with  H.  Montes,  G.  Rao  and  E.  Willis1  J 

(91 

and  with  H.  Montes  of  Isotopes,  Inc.  These  are  studies  of 

(8) 

observed  arrivals  from  Apollo  rocket  exhausts^  *  (long  period 
acoustic)  and  of  naturally  occurring  background  gravity  wave 


uclivily  ,  bulb  carriml  util  on  liio  1  but  opus  Doppler  lonobondu 
systom.  These  studies  indicate  that  propagation  and  background 
data  at  ionospheric  heights  can  help  refine  our  knowledge  of 
upper  atmospheric  structure  (neutral  gas  component) . 
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Abstracts  of  papers  in  press: 

1.  Phase  Height  Fluctuations  in  the  Ionosphere  between  130  and  250  km. 

I.  Tolstoy  and  II.  Montes 

Journal  of  Terrestrial  and  Atmospheric  Physics,  1971. 

Power  spectra  of  CW  phase-path  sounders  operating  at  6.0,  4.8  and 
2.4  MHz  often  display  peak  or  slope  changes  at  periods  in  the  5-20 
minute  range.  Plots  of  these  periods  vs.  real  height  tend  to  follow  mean 
VUisdia  frequency  curves  calculated  from  Standard  Atmosphere  models. 

This  suggests  that  statistical  properties  of  CW  phase-path  sounder  fluc¬ 
tuation  can  be  used  to  monitor  certain  mean  physical  parameters  of  the 
upper  atmosphere. 

2.  Generation  of  long  Internal  Gravity  Waves  in  Waveguides  by  Rising 

Buoyant  Air  Masses  and  other  Sources. 

I.  Tolstoy  and  J.  Lau 
Journal  of  the  Atmospheric  Sciences,  1971. 

The  displacement  fields  generated  in  an  internal  gravity  wave  waveguide 
between  plane  rigid  walls  are  compared  for  two  types  of  source:  an 
explosive  point  source  and  a  rising  buoyant  sphere  moving  at  constant 
speed.  It  is  concluded  that,  for  large  enough  spheres  and  comparable 
energy  expenditures,  the  buoyant  sphere  is  a  far  more  efficient  source 
of  long  internal  gravity  waves.  In  particular  it  appears  possible  to 
conclude  that,  in  the  case  of  large  events  such  as  nuclear  or  volcanic 
explosions  in  the  atmosphere,  the  rising  heated  air  mass  can  generate 
long  wavelength  (A  >  500  km)  internal  gravity  waves  at  ionospheric 
heights. 


3.  Dumping  of  very  long  Internal  Gravity  Waves  in  the  Atmosphere. 

I.  Tolstoy 

Journal  of  the  Atmospheric  Sciences,  1971. 

Assuming  most  of  the  energy  absorption  to  take  place  in  a  viscous  surface 
layer,  it  is  possible  to  make  a  simple  calculation  for  the  attenuation 
of  long  internal  gravity  waves  in  the  atmospheric  waveguide.  This  cal¬ 
culation  indicates  that  wavelengths  longer  than  a  few  hundred  kilometers 
may  propagate  to  considerable  ranges  as  waveguide  modes. 
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GENERATION  OF  LONG  INTERNAL  GRAVITY  WAVES  IN  WAVEGUIDES 


BY  RISiNG  BUOYANT  AIR  MASSES  AND  OTHER  SOURCES 


by 

I.  Tolstoy  and  J.  Lau 


Abstract 

The  displacement  fields  generated  in  an  internal  gravity  wave  wave¬ 
guide  between  plane  rigid  walls  are  compared  for  two  types  of  source:  an 
explosive  point,  source  and  a  rising  buoyant  sphere  moving  at  constant  speed. 
It  is  concluded  that,  for  large  enough  spheres  and  comparable  energy  expendi¬ 
tures,  the  buoyant  sphere  is  a  far  more  efficient  source  of  long  internal 
gravity  waves.  In  particular  it  appears  possible  to  conclude  that,  in  the 
case  of  large  events  such  as  nuclear  or  volcanic  explosions  in  the  atmos¬ 
phere,  the  rising  heated  air  mass  can  generate  long  wavelength  (X  >  500  km) 
internal  gravity  waves  at  ionospheric  heights. 


I.  Introduction 


Various  types  of  internal  gravity  wave  Green's  functions  have  been 
studied  in  the  geophysical  and  hydrodynamic  literatures.  Explosive 
sources  have  been  at  a  premium  in  that  particular  body  of  writing  which 
concerns  itself  with  the  atmospheric  effects  of  volcanic  explosions 
(Pekeris,  1948;  Scorer,  1950)  and  nuclear  tests  (Harkrider,  1964;  Hark- 
rider  and  Wells,  1968;  Pierce,  1965  and  1968).  That  type  of  investigation 
has  concerned  itself  primarily  with  acoustic  and  internal  gravity  wave 
modes  with  periods  less  than  ten  minutes  or  so.  Several  discussions  of 
the  effects  of  vertically  moving  sources  have  been  published  during  the 
last  decade  (Mowbray  and  Rarity,  1967;  Lighthill,  1967).  Of  particular 
interest  is  a  study  by  Warren  (1960),  giving  the  wave  resistance  for  an 
upward  moving  body  in  an  infinite,  incompressible  half-space;  his  results 
can  be  used  for  estimating  the  wave  amplitudes  generated  by  rising  hot 
air  masses.  The  possible  importance  of  this  type  of  mechanism  in  con¬ 
nection  with  nuclear  explosions  has  been  suggested  by  Knabe  and  Kahalas 
(1968). 

In  the  analysis  given  below  we  have  confined  ourselves  to  the  highly 
idealized  problem  of  a  perfect  internal  gravity  wave  waveguide  between  rigid 
walls.  This  problem  is  clearly  of  some  interest  to  oceanographers  for  it 
has  long  been  recognized  that  a  free  surface  is  approximately  s.  node  for 
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internal  gravity  waves  (.Lamb,  1945;  liekurt,  I960;  Tolstoy,  1963,,.  The 
assumption  of  a  rigid  top  surface  eliminates  the  surface  gravity  mode  while 
modifying  the  internal  modes  only  slightly,  whereas  the  en:.uing  analytical 
simplifications  are  considerable  and  make  this  approximation  worth-while. 

A  study  of  this  nature  is  germane  also,  to  some  extent,  to  a  class- of 
observed  atmospheric  disturbances  that  have  been  plausibly  ascribed  to  long 
wavelength  internal  gravity  wave  effects.  That  this  is  so  can  be  seen  as 
fol lows . 

The  upper  atmosphere  grades  continuously  into  the  "vacuum"  of  inter¬ 
planetary  space.  Well  below  the  height  at  which  the  mean  free  path  £.  of  the 
air  molecules  is  of  the  order  of  the  wavelength  X  of  a  hydrodynamic  wave, 
the  atmosphere  is  essentially  a  continuous  fluid.  Above  this  height,  i.e., 
for  X  <<  l ,  we  have  for  most  practical  purposes  a  vacuum.  Between  this 
vacuum  and  the  continuous  fluid  below  there  is  a  transition  region  in  which 
A  »  1.  To  understand  the  behavior  of  this  transition  region  one  should, 
strictly  speaking,  use  a  molecular  model  with  Boltzmann's  equation.  From 
the  macroscopic  point  of  view  this  region  is  that  in  which  the  kinematic 
viscosity  becomes  extremely  large:  it  is  this  layer  that  effectively  absorbs 
the  energy  of  upward  traveling  waves.  We  may  assume  (Tolstoy,  1967)  that 
the  thickness  d  of  this  transition  is  of  the  order  of  a  scale  height  and  is 
independent  of  A,  But  then,  for  A  >>  d  it  seems  reasonable  to  assume  that 
the  atmosphere  is  effectively  topped  by  a  free  surface  with  a  thin  (compared 
to  A),  high  viscosity  layer  just  under  it,  nnd  it  is  possible  to  imagine  long 
wavelength  waves  reflected  from  this  surface  with  modest  amounts  of  atten¬ 
uation.  Long  surface  gravity  waves  may  also  be  connected  with  this  'feffective 
free  surface"  (Tolstoy,  1967;  Tolstoy  and  Pan,  1970).  Furthermore,  as  has 
been  demonstrated  by  Harkrider  and  Wells  (1968),  the  low  order  long  wave¬ 
length  gravity  modes  are  often  channelled  at  heights  between  100  and  300  km.; 


3 


thus  perhaps  only  a  fraction  of  the  upward  travelling  energy  actually 
reaches  the  viscous  transition  layor  near  the  effective  "top"  of  the  atmos¬ 
phere  and  reflection  '..ill  be  close  to  total.  Observation  of  occasional 
round-the-world  atmospheric  gravity  wave  paths  suggests  that  reflexions  do 
take  place  in  the  upper  atmosphere:  certainly,  propagation  to  distances  of 
104  to  5  x  104  km.  (Breitling  ct  ol.,  1967;  Rose  et  al.,  1961;  Hultquist  et 
al.,  1961;  Dieminger  et  al.,  1962;  Herron  and  Montes,  1970;  Tolstoy  and 
Herron,  1970)  suggests  that  waveguide  effects  'take  place  for  these  long 
wavelengths,  just  as  they  have  been  demonstrated  to  occur  at  lower  heights 
for  the  shorter  period  acoustic-gr.  /ity  modes  (Harkridcr,  1964).  If  then  we 
confine  our  interest  to  the  internal  gravity  modes,  we  may  replace  the 
"effective  free  surface"  by  an  effective  rigid  surface  without  incurring 
serious  error  (see  comparative  group  velocity  calculations  for  a  free  and  a 
rigid  top  surface  on  a  layered  atmosphere  by  Tolstoy  and  Pan,  1970).  The 
propagation  of  long  wavelength  (A  >>  d)  internal  waves  in  the  atmosphere  as 
well  as  in  the  ocean  may  therefore  be  in  some  ways  similar  to  that  in  a 
waveguide  with  rigid  walls.  To  be  sure,  this  can  in  no  way  be  offered  as  a 
model;  however,  a  perfect  waveguide  of  this  kind  does  provide  a  simple  analog 
in  the  sense  that  many  of  the  long  range  propagation  effects  must  be  quali¬ 
tatively  similar. 

Our  somewhat  academic  looking  problem  should  thus  provide  us  with  clues 
concerning  the  propagation  to  great  ranges  of  long  internal  gravity  waves  in 
the  atmosphere  and  oceans.  We  are  going  to  show  that,  in  particular,  it 
gives  an  estimate  of  the  relative  efficiency  of  various  kinds  of  source, 
using  simple  analytical  approximations .  Our  analysis  suggests  a  possible 
explanation  for  the  efficient  generation  by  atmospheric  nuclear  tests  of 
long  internal  gravity  waves  in  the  ionosphere.  Such  waves  have  been  reported 
repeatedly  in  the  geophysical  literature  in  connection  with  ionosonde  obser- 
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vat  ions  of  large  yield  tests  and  arc  known  to  have  propaguted  to  very  great 
distances;  whereas  explosive  "pressure"  or  "energy"  sources  can  hardly 
account  for  the  magnitude  of  these  disturbances,  it  appears  possible  to  do 
so  with  u  "buoyant  rise"  mechanism  involving  a  rapidly  rising  mass  of  hot 
air. 


2.  Fundamental  equations 

In  the  following  we  shall  apply  the  method  of  normal  coordinates  to 
derive  the  wave  fields  originating  from  various  types  of  source.  This 
method,  used  extensively  in  acoustics  by  Rayleigh  (1945)  and,  later,  in  e.m 
field  theory  (Heitlcr,  1954),  can  be  applied  to  conservative  mechanical  or 
electro-mechanical  continua  of  any  sort.  As  pointed  out  by  Biot  and  Tols¬ 
toy  (1957)  this  method  has  several  advantages,  not  the  least  of  which  is  to 
provide  automatically  the  correct  orthonormality  relations  for  any  type  of 
wave  problem.  These  relations  (equation  17)  occur  naturally  in  terms  of 
the  displacement  eigenfunctions  and,  as  a  result,  the  formalism  using  dis¬ 
placement  fields  is  particularly  suitable.  Thus  we  begin  with  the  free 
field  equations  in  the  form  (Tolstoy,  1963): 

pg  -  VAe  +cgV;  -  lzPgE  ■  0  (1) 

where  g  is  the  perturbation  displacement  of  the  stratxf 1H  fluid,  having  r, 
z  components  £,t;  is  the  unit  vector  in  the  z  direction,  p  is  the  unper¬ 
turbed  equilibrium  density,  g  the  acceleration  of  gravity,  and 

(2) 


c  ■  v*d  *  if  +  UL  +  iLl 
*  r^  3r  3z 

and  A  is  Lame's  constant 


A  =  pc^  (3) 

c  being  the  speed  of  sound.  In  these  equations  it  is  to  be  assumed  that  A,p 


are  functions  of  z  only. 

Assuming  a  constant  c  with  an  exponential  density  p: 


S 

-2vz 

(4) 

P  ■  P0e 

The  solutions  of  (1)  have  the  form: 

;  ■  evzJo(<r)h 

(5) 

(6) 

with 

f  .  ,cb“2  ~  +  (v  -  i&)h 

i 

(7) 

2 

u2  w  2 

»  — 7-  - 

cz 

.  (8) 

h  obeys  the  Helmholtz  equation: 

d2h  2u  n 

'  dP-  *  r  h  =  ° 

(9) 

Y2*  b2  +  -^jN2  -  v2 
ur 

(10) 

Here  N  is  the  Vaisala  frequency,  i.e.. 

N2  =  2vg  -  £ 

(11) 

In  this  case  then, y2  is  not  a  function  of  z. 

The  modes  of  a  waveguide  of  thickness  H,  with  rigid 

walls  at  z  ■  0, 

z  ■  H  will  then  be: 

c  -  vvz  sinv  V<r> 

(12) 

«  "  VV*  F*  IYmC0SV  +  (V  "  cfe5inYB*]Ji(fr) 

(13) 

with 

y^H  ■  mu  ,  ra  »  0,  1,  2,  .  .  . 

(14) 

Thus  (12),  (13)  are  of  the  form: 

d  ■  q_a 
=m  Tn'm 

(15) 

q_  is  the  normal  coordinate  ard  is  a  function  of  t  only. 

obeying  the 

oscillator  equation.  In  the  absence  of  exciting  forces: 

0 


(16) 


q  ♦  urq  = 
nn  in'm 

’•  The  orthonormality  conditions  (Biot  and  Tolstoy,  1957)  are: 


Pa  *a  d t  =  6 
vm  *n  mn  m 


(17) 


where  the  integration  extends  over  the  whole  of  physical  space  and  is 
the  Kronecker  delta: 


6  =1  m  =  n 

om 


(18) 


»  0  m  ^  n 

In  applying  (17)  to  (12),  (13)  and  infinite  spaces  one  must  use  symbolic 
results  (Biot  and  Tolstoy,  1957;  Tolstoy  and  Clay,  1966)  of  the  type: 


*>  roa 

J0(<r)rdr  -  J*0cr)rdr  ■ 

'  n 


(19) 


'c  'o 

Thus  the  solutions  (12),  (13)  give: 
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or 


with 
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°m  =  1  +  +  (v  “ 


(21) 


(22) 


In  the  presence  of  exciting  forces,  the  will  be  the  solutions  of: 

,2 „  =  Sm  =  -  KdK.Qm 
>oH  °m 


%  +  ^  "  wpjT’On 


(23) 


Here  is  the  component  of  generalized  force,  deduced  from  the  actual 

force  F  by  the  principle  of  virtual  work  (Biot  and  Tolstoy,  1957) : 


F*a  di 
*  *m 


(24) 
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If  then  we  say  that  G  (t)  is  the  solution  of 

m 

G  ♦  ui2G  -  (25) 

will  have  the  forms: 
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(26) 


£  3  jz'1  u'  eVZ  l  Ty  cosy  z  +  (v  — $-)siny  zl 
up  H  L  t'm  v  'm  J 

o  m 


G 


A--J2J.(<r)KdK(27) 
°  m 


Note  that  the  solution  of  (25)  is,  in  most  cases  of  interest,  available  in 
closed  form  after  applying  the  convolution  theorem: 


m 


€ 


Q(T)sinw(t  -  T)dr 


(28) 


Thus  the  t,  £  field  will  be  available  through  (26),  (27)  by  means  of  a 
single  integration  which  can  usually  be  carried  out  with  the  help  of  the 
method  of  stationary  phase  (Lighthill,  1965). 

In  this  paper  we  shall  compare  two  types  of  source: 

1.  An  explosive  source,  corresponding  to  an  instantaneous  injection  of  volume 

2.  An  upward  moving  force,  providing  a  crude  simulation  of  a  rising  hot  air 
mass  of  dimensions  small  compared  to  the  wavelengths  being  studied. 

A  simple  representation  of  the  explosive  source  is  given  by  a  force 
.compressing  a  small  sphere  which  is  suddenly  released  at  t  ■  0.  In  other 
words  we  have  a  force  of  time  dependence  1  (- 1)  and  we  shall  write: 

^  <29) 

i.e. ,  by  (28)  : 

Gm  s  •  C0Sut)Ajn  (30) 

In  the  problem  under  discussion  one  is  interested  only  in  the  time  dependent 
terms,  and  we  may  take: 


b 


G  = 
m 


COSuit  . 

77- A 


(31) 


or  m 

It  can  be  shown  (Biot  and  Tolstoy,  1957;  Tolstoy  and  Clay,  1966)  that  the 
"unit  explosion"  corresponding  to  the  injection  of  a  unit  volume  at  t  ■  0 
gives: 


A  =  -X7*am  (32) 
m  =>m 

where  all  quantities  are  to  be  evaluated  at  the  source  point  r  =  0,  z  *  zQ. 
If  then  B  is  a  measure  of  the  source  strength  (B  =  1  corresponds  to  a  unit 
volume)  we  write: 


Qm  "  IV2cosyiZo  *  (vu)2  ‘  *28)sinV0]l(-t)  (33) 

And  by  (29)  -  (31)  we  shall  have  in  (26),  (27): 

Gm  *  Bpoe"V“#p'  tV2c0SVo  *  Cv“2  ‘  (34) 

We  thus  have,  in  integral  form,  an  explicit  expression  for  the  field  due 
to  an  explosive  point  source  at  r  »  0,  z  «  zq. 

The  other  type  of  source  we  are  interested  in  is  a  concentrated 
force  switched  on  at  t  «  0,  z  ■  zq,  moving  upward  with  the  steady  velocity 
V  and  switched  off  at  t  ■  T  so  that  the  total  distance  travelled  is 

L  *  VT  <  H  (35) 

For  a  vertical  force  of  this  kind, 

E  *  lz6(r)6[(z  -  zQ)  -  VtJ(z)  (36) 


We  may  envisage  two  somewhat  different  cases,  corresponding  to: 

F  -  Fj  =  F1 pH 


F  =  F2  =  F2p 


The  case  F  =  would  appear  to  correspond  to  the  long  wavelength  limit  of  a 
small  buoyant  sphere,  since  Warren  (1960)  has  shown  that  in  this  case  the 
wave  resistance  upon  the  rising  sphere  should  be  proportional  to  the  local 
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value  of  the  density.  In  fact,  there  is  reason  to  doubt  that  Warren's 
result  can  be  extended  to  the  violent  events  accompanying  an  atmospheric 
thermonuclear  explosion,  and  we  shall  show  below  that  at  this  stage  of  the 
art  at  least,  there  is  little  to  choose  between  the  two  models  (37),  (38). 
However,  we  shall  use  Warren's  analysis  to  deduce  an  order  of  magnitude  for 
F2,  F[  from  the  physical  parameters  of  a  hypothetical  rising  sphere  of  hot 
gas.  Note  that  F;  is  a  measure  of  energy  and  the  F2  source  corresponds  to 
an  energy  density  input  ii.to  the  wave  field  that  falls  off  with  height  like 

k 

p  ,  whereas  the  Fj  source  gives  a  constant  energy  input  with  height.  If,  in 
fact,  the  full  height  of  rise  L  does  not  exceed  p  scale  height,  the  two 
kinds  of  source  will  give  similar  orders  of  magnitude  insofar  as  the  field 
displacement  amplitudes  are  concerned. 

For  F  »  Fj  we  have,  by  (12),  (24),  (36): 

Q  ■  sin(u  t  ♦  <J>  )eVZ°  p*F.  0  <  *t  $  T 
in  om  m  o  i 


-  0 


T  <  t 


(39) 


where 


b)  =  Y  V 
om  m 


fr  ■  Y  z 
m  mo 


(40) 

(41) 


Thus  we  have,  driving  the  oscillator, a  harmonic  force  switched  on  at  t  »  0, 
switched  off  at  t  =  T.  The  convolution  theorem  gives  us  then,  for  t  >T: 

„  eVI°(&F1  fT 


m 


u 


sin(w  t  +  $  )sinw(t  -  r)di 
om  m 


(42) 


i.e. 


G  *  F  —e  $  {— - sir^r(b)  +  u^m)cos[bit  -  <p  -  4(u  ♦  u _ )] 

m  lb)  o  u+b) .  _  2  om  L  m2  om  J 

om 


(43) 


IT  T 

-  sin^-(b)  -  bi)  COS  [bit  ♦  4  +  Tr(b)  -  bl)l) 

bi  ~u  2  om  w  m2  om  J 
om 


On  the  other  hand,  we  have  for  F  ■  F2: 
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-vVt  .  ,  „  , 

e  sinfu  t  * 
om 


vz  n 
<p  )o  ° 
Ym 


*0F2 


0  <  t  s  T 


(44) 


=  0  T  <  t 

The  harmonic  force  driving  the  oscillator  in  the  time  interval  0  <  t  <  T 
is  now  modulated  by  the  factor  e"vVt.  The  corresponding  G,,,  is  easily 
derived,  but  we  shall  not  need  to  refer  to  it  and  will  not  write  it  out  here. 


3.  Approximations  and  the  displacement  field 

In  the  previous  section  we  have  obtained  the  displacement  wave  fields 
of  the  internal  gravity  and  acoustic  modes  generated  by  two  different  types 
of  source  in  an  isothermal  layer  of  fluid  between  plane  rigid  boundaries. 
These  solutions  (equations  26  and  27)  appear  as  a  series  summation  for  the 
displacement  field;  each  term  of  the  series  is  the  contribution  of  a  given 
waveguide  mode  and  has  the  form  of  an  integral  over  k.  To  put  these  results 
in  useful  form  one  must  evaluate  the  integral:  this  can  be  accomplished 
approximately  by  the  method  of  stationary  phase.  As  is  well  known (Lighthill, 
1965)  this  method  is  the  mathematical  embodiment  of  the  physical  fact  that 
the  energy  within  a  narrow  band  travels  with  the  group  velocity;  it  is 
equivalent  to  the  statement  that  the  principal  contribution  to  an  integral 
such  as  (26)  comes  from  the  vicinity  of  that  point  in  the  u,k  plane  for  which 

•~(wt  -  <r)  =  Ut  -  r  =  0  (45) 

where  U  is  the  r  component  of  the  group  velocity.  In  our  case  this  principle 
is  best  applied  to  the  "far  field"  solutions,  keeping  only  the  leading  term 
of  the  asymptotic  expansion  for  JQ(tc r) .  Thus,  if 

I  B  I  P(K)cos(wt  ♦  ip) JQ  (ter)  ted*  (46) 

'o 

we  shall  take 

*00 

I  =  2”‘*(irr)”M  P(»c)[cos(»cr  -  —  -  ut  -  t|»)  ♦ 

7T  h 

cos  (ter  -  ♦  ut  ♦  i^)]teAl»c 


(47) 
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Since  furthermore  U  represents  the  velocity  of  energy  transport  it  is  always 
positive,  i.e.,  away  from  the  source;  thus,  corresponding  to  (45),  we  keep 
only  the  first  term  in  brackets  in  (47)  and  write: 


I  =  2 


(Ttr)_isRe 


.00 


P(*)e 


i(<r  -  ut  - 


,/4)Au 


(4S) 


The  classic  result  of  the  method  of  stationary  phase  is: 

I  .  PU  )  1  V* - re'i(*  *  -  Oot  -  n/4sgnU ' )  (49) 

0  AV|% 

where 

U'  =  ||J  (501 

and  k0,uj0  are  given  by  (45)  as  explicit  functions  of  r,t. 

In  principle,  then,  we  are  in  a  position  to  evaluate  (26)  and  (27), 
with  Gm  given  by  (34)  (explosive  source)  or  (45) (buoyant  rise  source) ;  but 
to  utilize  the  result  (49)  we  must  first  solve  (45).  For  the  general  case 
of  a  compressible,  stratified  fluid  in  a  gravity  field  this  leads  to  an 
algebraic  equation  of  the  fourth  degree  in  u>2  or  <z  and  the  solution  is 

best  carried  out  numerically.  However  we  are  interested  mostly  jn 

\ 

establishing  orders  of  magnitude  and  in  some  of  the  general  properties  of 
the  internal  gravity  wave  modes,  so  that  wo  shall  use  an  additional  approxi¬ 
mation  which  will  allow  us  to  obtain  simple  analytical  results.  This  is 
the  assumption  of  incompressibility. 

The  validity  of  the  incompressible  approximation  hinges  upon  the  fact 
that,  in  an  isothermal  half  space  one  may  write,  for  the  internal  gravity 
wave  branch  and  for  given  (Tolstoy,  1963): 

u.2 

w2  3  w?(l  +  — ^2"  +  .  .  •)  (51) 

ua 

where  is  the  value  given  by  the  incompressible  formula  (while  keeping 
the  actual  N,v  values  characteristic  of  the  isothermal  model)  and  is  the 


value  of  ui  given  fox'  tho  acoustic  branch  after  neglecting  gravity.  For 
given  *, 


“a  va 


(52) 


where  v.,  v  are  the  corresponding  phase  velocities;  as  is  well  known, 

1  & 

-*■  oo  as  <  •*  0,  whereas  v^  remains  finite  and  for  long  wavelengths  the 
incompressible  approximation  will  be  quite  good  even  when  applied  to  gases. 
Thus  we  take: 


2  2[N2  , 

Y*  =  K*  — T  -  1 

UT 


-  V' 


(S3) 


where  N  has  the  numerical  value  corresponding  to  the  compressible  model 
(i.e.,  N  is  a  constant  given  by  equation  11).  We  are  thus  led  to  the 
results: 


U  »  U 


N(Vft  ♦  v2) 


m 


(Y 


2  +  y2  +  ^2^  fi 
m  ' 


(54) 


where  has  the  value  given  by  (14).  Substituting  (54)  in  (45)  yields  then 
for  each  mode: 


<  a  N  V  1  8* 
m  m 


(55) 


where  Vm  is  the  maximum  (low  frequency  limit)  of  the  group  and  phase  velocities 


of  each  mode: 

V  *  N(v2  +  y2)  ^ 
m  '  'm' 

and 

8  ■  t^3  r”2^3  V%  -  1 


m 


m 


likewise: 


and: 


ui  -  NV 
m  m 


-V3  t"  V3  j.  V3  qH 
m  m 


(56) 


(57) 


(58) 


u  t  •  k  r  »  rNV" 1 8  ^ 
mm  mm 


(59) 


Also  note  that  (22)  becomes,  by  virtue  of  the  incompressible  assumption 


N 


(60) 


m  " 

We  may  then  use  (49)  together  with  (55)  -  (60)  to  obtain  explicit,  analytic 
approximations  for  (26),  (27)  for  the  two  kinds  of  source  corresponding  to 
(34)  and  (43). 

Thus,  for  the  internal  gravity  modes  generated  by  an  explosive 
source  we  have  the  vertical  displacement  field: 

ir"1H"13"JSN-1eV(z'z°^B  l  [gV'^t^V^sinY  z 

m  mo 

m 

(61) 

"  Vm/3t"  /3r"  /3(ymCOSVo  +  VsinYmzo}  -nY-z  cosrNV^e/2 
Whereas  the  rising  force  Fj  gives: 

C  -  Tr1ir13"3sp-1r_1eVCz+zo)F  £  V" 1  e^: 
buoy  o  1  u  m  m 


’expl 


m 


m  m 


m 


sinY  z 
mm  m 


(62) 


.  T 
smxy 


L  x 


cos  u  - 


.  T 
smyy 

y 


cos  v 


with 


-  u  ♦  NV‘1/3eV1/3r^3 


om  mm 


u  -  Nv-1/%\‘^3r^3 


om  m  m 


u  »  xr  +  Y  z  -  rNV_16 
2  mo  mm 

v  =  yJ  ♦  Y  z  +  rNV_10^2 
72  mo  mm 


(63) 

(64) 

(65) 

(66) 


In  these  results,  each  mode  starts  at  t  *  rV'1:  the  amplitude  is  zero  for 

times  t  <  rV"1. 

m 

Note  that  for  constant  r,  as  t  a  |c|  -*■  09  like  t  ^3.  In  practice, 


however,  this  need  not  disturb  us  for  it  is  the  same  kind  of  divergence 
that  occurs  in  the  analogous  Cauchy-Poissun  problem  for  surface  gravity 


H 


waves  (Lamb,  1945) ;  it  is  due  to  the  concentrated  6 (r)  nature  of  the 
exciting  force  and  may  be  eliminated  by  taking  a  source  of  finite  width. 

We  are  interested  here  in  the  long  wavelength  early  arrivals,  for  U  *  300  m. 
sec"1  or  so,  for  which  the  wavelengths  are  much  greater  than  the  width  of 
the  actual  source,  so  that  we  may  ignore  this  effect. 

The  results  (61)  and  (62)  correspond  to  the  highly  idealized  case  in 
which  the  attenuation  of  the  waves  is  neglected.  In  practice,  internal 
gravity  waves  in  planetary  oceans  and  atmospheres  will  be  attenuated  by  a 
variety  of  mechanisms.  If  we  wish  to  make  our  results  in  any  sense  analo¬ 
gous  to  cases  liable  to  occur  in  practice,  we  must  introduce  attenuation  in 
some  form.  In  fact,  it  is  generally  believed  that  the  chief  mechanism  of 
attenuation  of  loi.g  wavelengths  in  the  atmosphere  is  due  to  viscosity  and, 
in  particular,  to  the  highly  viscous  "transition"  layer  mentioned  in  the 
introduction;  i.e.,  near  the  effective  "top"  of  the  atmosphere  the  attenu¬ 
ation  increases  to  such  large  values  that  all  short  and  moderate  wavelengths 
reaching  these  heights  are  completely  absorbed.  However,  as  pointed  out  in 
the  introduction,  large  atmospheric  explosions  have  been  known  to  generate 
long  wavelength  disturbances  corresponding  to  displacements  of  several  km. 
at  ionospheric  heights  and  traveling  out  to  distances  of  10*4  to  5  x  10*4  km. 
(i.e.,  around -the- world  paths).  These  long  paths  suggest  that  on  such 
occasions  gravity  waves  are  ducted  and  must  in  some  manner  be  reflected, 
either  at  an  "effective  free  surface"  or  both  there  and,  partially,  at  a 
lower  level;  whatever  the  precise  mechanism  of  the  ducting,  the  experimental 
results  suggest  that,  on  occasion,  the  attenuation  in  the  horizontal 
direction  is  small  enough  to  allow  for  a  number  of  reflections  before  the 
wave  dissipates  too  much  of  its  energy.  Thus  we  may  apply  the  observation 
(Lighthill,  1965)  that  small  amounts  of  attenuation,  although  important  in 
practice,  do  not  change  the  group  velocity  and  stationary  phase  theorems;  it 
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is  therefore  possible  to  take  the  attenuation  into  account  by  simply  multi- 

-  <5mr 

plying  the  undamped  solutions  by  e  m  ,  where  6m  is  some  function  of  the 
frequency  or  wuvenumber  dictated  by  the  physics  of  the  process.  Since  we 
are  dealing  here  with  a  viscous  mechanism,  we  may  assume  that  6^  is  propor¬ 
tional  to  u2  (at  least  for  order  of  magnitude  calculations,  when  the  phase 
velocity  does  not  vary  widely  over  the  frequency  band  of  interest).  Thus, 
for  added  realism  in  our  solutions  (61),  (62),  we  may  multiply  each  mode  by 

e-nu.2r  s  e-nV^3t'2/3r2/3era  (6?) 

Sensible  numerical  values  of  the  constant  n  may  be  secured  on  the  assumption 
that  the  wavetrain  will  give  detectable  amplitudes  at  some  range:  like,  say, 
an  order  of  magnitude  decrease  for  10**  or  5  x  104  km. 

4.  Results  and  discussion 

Comparing  the  -solutions  (61)  and  (62)  allows  us  to  demonstrate  the  fact 
that  a  buoyant  rise  type  of  source  (such  as  Fj  or  F2)  is  much  more  effective 
in  exciting  internal  gravity  waves  than  an  explosive  source,  at  least  in 
connection  with  large  energy  sources  such  as  volcanic  explosions  and  thermo¬ 
nuclear  tests  which  create  a  big  enough  volume  of  rising  hot  fluid. 

In  view  of  Warren's  results  (1960)  showing  that  the  wave  resistance  on 
»  buoyantly  rising  sphere  is  proportional  to  the  density  of  the  surrounding 
medium,  it  would  be  logical  to  assume  that  it  is  the  F2  source  that  simulates 
a  rising  sphere  of  radius  small  compared  to  the  wavelengths  of  interest. 

In  fact  though,  the  generalized  forces  (39)  and  (44)  will  give  results  of 
the  same  size  for  moderate  rise  speeds  and  times  V,T,  i.e.,  as  long  as  vVT 
is  not  too  large.  Certainly,  if  we  were  to  take 

Fl  •  F2<pl fi>  (68) 

where  the  average  is  taken  over  the  rise  height  L,  we  may  expect  the  Fj  and 
F2  sources  to  give  similar  orders  of  magnitude  for  the  displacement  ampli- 
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tudcs  (even  though  t lie  details  of  t he  field  may  bo  different). 

Adopting  (68)  then,  wo  can  estimate  F2  by  assuming  that  1’2  is  essen¬ 
tially  the  total  wave  resistance  to  the  upward  motion  of  the  sphere  and 
write : 

F2  =  Rd  x  2vga4  .  (69) 


where  R^  has  been  given  6y  Warren  as  a  function  of  the  radius  a  of  the 
rising  sphere.  For  radii  that  are  not  excessive  (Warren,  1960): 

<*  a22vgn_2V“2  (70) 

so  that 


F2  *  4ir‘2V"2v2g2a6  (71) 

The  effective  radius  a  of  the  sphere  depends  upon  the  mechanism  involved. 

In  the  case  of  thermonuclear  atmospheric  explosions,  this  may  be  taken  as 
(Pierce,  1968) 

a  =  0.4  eVV3  (72) 

where  E  is  the  total  energy  of  the  explosion  and  p  is  the  atmospheric 
pressure  at  z  «  ifl.  Using  this  result  togother  with  (71)  gives 

F2  »  1.6  x  10-3v2g2p-2V”2E2  •  .  (73) 

S\ 

V  also  depends  upon  the  energy,  but  in  a  manner  hard  to  ascertain.  If  we 
use  the  buoyant  rise  velocity  formula  given  by  Scorer  (1957)  for  less  extreme 
conditions,  we  take: 


V  =  (gfla)*5 

where  8  is  the  buoyancy  defined  as 


(74) 


(75) 


Ap  being  the  difference  between  the  densities  outside  and  inside  the  rising 
sphere.  It  is  very  difficult  to  estimate  8  on  theoretical  grounds,  but  if 
we  assume  that  for  the  very  great  temperatures  involved  in  a  thermonuclear 
explosion  8  9  1,  then 


E.g.,  if  zq  «  35  km.,  p  «  102  kg.m"2  and  if  we  assume  that  for  1  megaton 
E  a  1013  kg.m,  we  have  V  «  102  m  sec-1. 

Using  (76)  we  rewrite  (73)  as 


Fa  4  x  10"2v2  E5/3p_5/3 


(77) 


where  we  have  assumed  g  =  10  m.sec"2.  This  result  suggests  that  F2  is 
actually  an  increasing  function  of  zq.  However  we  know  that,  in  practice, 
the  hot  sphere  slows  down  and  reaches  a  ceiling  and  then  starts  spreading: 
thus  V  -*•  0  ar.d  at  this  point,  according  to  (73),  F2  would  increase  indefi¬ 
nitely  which  is  clearly  impossible..  In  fact  F2  must  decrease  to  0  as  the 
air  mass  comes  to  a  halt.  In  view  of  these  complicating  factors  it  appears 
reasonable  to  use  a  constant  energy  input  source  of  the  type  :  pending 
a  more  quantitative  analysis  of  the  problem  we  assume  that  this  will  at 
least  give  us  tiie  proper  orders  of  magnitude. 

For  order  of  magnitude  calculations,  then,  we  assume  an  F}  source,  with 
Fj  given  by  (68).  The  quantity  F2  in  this  equation  is  then  calculated  from 
(77)  using,  for  example,  the  actual  average  atmospheric  parameters  of  the 
lower  100  km.  of  the  atmosphere  (v  »  7  x  10" 5)  rather  than  those  of  the 
waveguide  model  (it  is  the  magnitude  of  the  source  we  are  modelling  at  this 
point).  In  the  analog  waveguide  it  is  probably  more  consistent  in  (68)  to 
use  the  v,p  values  of  the  waveguide,  i.e., 

vVfl 


Fi  *  F2 


1-e 


vVT 


-vzf 


(78) 


where  V  is  an  assumed  rise  rate  compatible  with  (76),  T  the  duration  of  the 
rise,  v  is  characteristic  of  the  model  (smaller  than  that  of  the  lower  100 
km.  of  the  atmosphere  and  larger  than  that  of  the  upper  atmosphere).  Thus, 
for  V  ■  100  m.sec"1,  T  ■  300  sec.,  v  “  3  x  10"5  we  have,  for  order  of 
magnitude  purposes: 
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F i  =  F.,  o'vZ°  x  0.  7  (79) 

Insofar  as  the  explosive  source  is  concerned,  we  may  recall  that,  for  a  • 
unit  injection  of  volume,  B  *  1  (Biot  and  Tolstoy,  1957).  This  corresponds, 
approximately,  to  1  kg.  of  conventional  chemical  explosive  which,  upon 
detonating  at  sea  level  atmospheric  pressure,  is  transformed  into  1  m3  of 
gas  and  gives  an  energy  release  of  the  order  of  104  kg.  m.  At  higher 
altitudes,  the  volume  injected  in  this  manner  by  a  given  charge  size  grows 
like  p*1,  i,e.,  like  e2vZ°.  On  the  other  hand,  the  effective  volume 
injected  in  this  model  does  not  grow  proportionately  to  the  mass  of  the 
explosive:  the  detonation  wave  travels  through  the  explosive  and  the 
effective  injected  volume  is  much  smaller.  Indeed,  the  usual  scaling  for 
amplitude  -  and  thus  for  B  *  is  the  cube  root  law.  Thus: 

B  =  E1/a  10_4/3  e2vZ°  (80) 


Using  (80),  (79)  with  (61),  (62)  we  may  determine  the  relative  efficiency 
of  the  buoyant  and  explosive  sources.  We  are  interested  in  the  earlier 
arrivals  for  which  Um  »  rt"l  =  Vm;  so  that  we  may  write,  approximately,  for 
g  ■  10*  N  «  1.5  x  10"2: 


|;|  1  *  (r-^-^-^'Sr-1  ev(z*Zo))  x  2.6  x  10  x  E^3  '  (81) 

Note  that  the  explosive  arrival  starts  with  a  discontinuous  step  at  t  ■  rV 
for  each  mode. 

The  buoyant  arrival,  on  the  other  hand,  is  seen  from  (62)  to  rise 
gradually  from  zero  amplitude  at  t  =  rV"1 .  However,  it  is  easily  verified 
that  fully  developed  amplitudes  correspond  to  0^  i  10" 1  (see  figures 
for  example).  Thus: 


.1 

m 


or, 


k.  I  2  {r-17T-1H-13-isV- 1 
'*buoy'  m 


for  T  »  3  x  102  sec: 

|c|,  2  (r“lir*1H-13JsV-1 

1  1  buoy  m 


I  e*8  e'VI°  x  0.7  x  F, 

t  m  * 

ev(,+,°)>  7.4  x  10-17  x  e7/3VZOEV3 


(82) 


(83) 


iii 


ill  us,  the  ratio  of  the  maximum  amplitudes  for  the  early,  long  wavelength, 
internal  gravity  wave  trains  excited  by  the  buoyant  rise  and  explosive 
mechanisms  is,  at  least 


R 


*buoy 


*expl 


-  Z  x  10-18 


a) 

p  . 


7l 


>/3 


(84) 


where  E  is  the  equivalent  energy  in  kgm.  Thus  the  relative  efficiency  of 

the  buoyant  rise  mechanism  increases  rapidly  with  the  energy  (E  and 

_  7/ 

with  the  height  of  initiation  (p  /6) .  The  estimate  (84)  shows,  for 
instance,  that  for  E  =  10 12  kgm  (100  kT) ,  and  p  =  10"3pQ  (z0  =  km) , 

R  z  102.  Thus  large  explosions  at  moderate  to  great  heights  have  a  ten¬ 
dency,  through  the  buoyant  rise  mechanism,  to  be  efficient  generators 
of  long  wavelength  ionospheric  gravity  waves:  the  explosive  phase  of 
the  mechanism,  on  the  other  hand,  is  a  relatively  negligible  source  of 
such  waves.  This  should  hold  equally  well  for  volcanic  explosions. 

Insofar  as  the  absolute  amplitudes  of  excitation  of  the  internal 
gravity  wave  modes  by  the  buoyant  rise  is  concerned,  it  is  unlikely 
that  the  above  formulae  would  give  results  pertinent  to  the  real  atmos¬ 
phere.  However,  one  may  point  out  that,  at  least,  these  results  suggest 
that  very  large  amplitudes  can,  in  principle,  be  obtained.  For  instance, 
figure  1  shows  that  for  r  -  104  km,  assuming  H  2  5  x  105  m.,  and  for  a 
1  MT  explosion  at  zQ  «  100  km,  vertical  amplitudes  in  excess  of  10  km 
are  easily  justified.  In  fact,  our  model  underemphasizes  the  displace¬ 
ment  amplitudes  to  be  obtained  at  great  heights  because  it  uses  an  average 
exponential  density  law;,  in  addition  we  know  that  ducting  effects  occur 
for  these  wavelengths  above  100  km. 

Figure  1  shows  a  succession  of  waveforms  for  c  calculated  for  the 


model: 


II  =  5  x  J.O5  m 
N  =  1.5  x  10-2  rad.  sec"1 
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v  =  3  x  10" 5  m- 1 
V  =  102  m.  sec-1 
T  =  3  x  102  sec. 
at 

r  =  107  m. 
z  =  2.5  x  105  m. 

and  for  several  zQ  values  and  various  attenuation  coefficients  n 


(85) 
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We' hove  c-stabli  shecf  a  suitable  Green's  function  model  to  explain  the  generation  of 
lung  period,  large  amplitude  gravity  modes  by  thcrmorcic lenr  and  nuclear  tests  in  the 
atmosphere.  Thus,  it  is  shown  in  a  paper  in  press  by  Tolstoy  and  Lau,  that,  the 
important,  mechanism  here  must  be  the  rising  shore  of  hot  gas .  It  is  not  difficult 
to  demonstrate  that  ••  for  large  yields,  at  least  -  this  mechanism  should  be  orders 
of  t:v,\  u  i  tede  more  efficient  in  generating  long  period  internal  gravity  wave  modes 
efficient  in  gcncratip^  long  period  internal  gravity  wave  modes  than  the  explosion 
phase  proper.  It  is  our  belief  now  that  this  is  the  source  of  observed  long  wave¬ 
length  gravity  waves;  although  wc  have,  in  the  attached  report,  merely  illustrated 
the  principle  of  the  mechar i sm ,  without  extending  the  results  to  a  realistic  model 
of  the  atmosphere ,  tiu-.ro  appears  to  be  little  doubt  that  one  may  in  this  manner  . 
account  for  the  large  r.n.-pl  i  t.udes  of  vertical  motion  observed  for  such  disturbances 
traveling  in  the  j  one  spilt- r  >.  . 
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Figure  Caption 


Figure  1 

Internal  gravity  wave  vertical  displacement  amplitudes  (in  m)  at  height 
2  =  250  km  in  waveguide  model  of  equation  85,  excited  by  buoyant  hot  air 
mass,  rising  at  100  m  sec-1  for  five  minutes,  starting  from  various  heights 
z  »  zQ  (sum  of  first  seven  modes) .  Parameters  have  been  estimated  for  a 
1  MT  explosion.  Several  possible  attenuation  factors  h-ve  been  used  for 
qualitative  illustration  of  effects  of  different  attenuation  magnitudes. 
Although  these  calculations  can  in  no  way  be  construed  as  giving  a  model 
of  atmospheric  behavior,  the  results  do  suggest  that  this  mechanism  can 
create  large  amplitude  disturbances  at  ionospheric  heights. 


